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ABSTRACT 

We derive the distribution of neutral and ionized gas in high redshift clouds which 
are optically thick to hydrogen ionizing radiation, using published data on Lyman-limit 
and Damped Lyman-a absorption systems in the redshift range 1.75< z <3.25. We 
assume that the distribution of the hydrogen total (HI+HII) column density in the 
absorbers, Nh, follows a power law KNjj a , whereas the observed HI column density 
distribution deviates from a pure power law as a result of ionization from a background 
radiation field. We use an accurate radiative transfer code for computing the rapidly 
varying ratio Nh /Nhi as a function of Nh- Comparison of the models and observations 
gave excellent fits with Maximum Likelihood solutions for the exponent a and for X, 
the value of log(iV# /Nhi) when the Lyman-limit optical depth along the line of sight is 
Tjjjj = 1. The slope of the total gas column density distribution with its relative 3a errors 
is a = 2.7^ 7 an d X = 2-75 ± 0.35. This value of X is much lower than what would be 
obtained for a gaseous distribution in equilibrium under its own gravity. The ratio ryo of 
dark matter to gas density is however not well constrained since log(ryo) = 1.1 ±0.8. An 
extrapolation of our derived power law distribution towards systems of lower column 
density, the Lyman-a forest, tends to favour models with log(ryo) ~ 1.1 and a ~ 2.7 
to 3.3. With a appreciably larger than 2, Lyman-limit systems contain more gas than 
Damped Lyman-a systems and Lyman-a forest clouds even more. Estimates of the 
cosmological gas and dark matter density due to absorbers of different column density 
at z ~ 2.5 are also given. 
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Subject headings: QSO: absorbers — radiative transfer — dynamic: dark matter 

1. Introduction 

Lyman limit systems (hereafter LLS) are detected as intergalactic clouds which absorb the 
quasar radiation at energies above the ionization edge of neutral atomic hydrogen. The optical 
depth to the continuum radiation at the Lyman edge is tll = Nhi /(L6 x 10 17 cm -2 ), where Nhi 
is the HI column density along the line of sight. For tll ~ 1, Nhi can be measured accurately 
regardless of the velocity structure of multiple components; for tll 3> 1 the absorption can be 
detected, although Nhi cannot be measured. We shall not analyze in detail the Lyman-a forest 
systems, at much smaller Nhi, but will consider the whole range of column densities from the LLS 
to the, mostly neutral, damped Lyman-a systems (hereafter DLS) with much larger Nhi- We will 
then discuss a possible extrapolation of our results towards the Lyman-a forest region. There have 
been some suggestions that LLS, as well as Mgll absorbers with similar HI column densities, may 
be related to outskirts of galaxies (Bergeron & Boisse 1991; Gardner et al. 1999; Rao & Turnshek 
2000) while DLS are more likely to originate from their inner regions (Prochaska Sz Wolfe 1998, 
and references therein). We shall show that, whether this suggestion is verified or not, the ratio 
of the total to neutral hydrogen column density for typical LLS can be estimated assuming only 
that the total gas column density distribution g(Nn±) is a smoothly varying function from LLS to 
DLS range. In particular we shall derive a value for X, the logarithm of Nh±/Nhi±_ for a line of 
sight value of Nhi = 1-6 x 10 17 cm~ 2 , at intermediate redshifts, namely 1.75 ^ z ^ 3.25. From X 
we shall derive a, the exponent in the power law for g. 

If the dimensionless parameter X is very large, the ionization fraction changes rapidly and the 
value of Nhi can increase from ~ 10 17 cm -2 to ~ 10 19 cm~ 2 over a fairly small increase of the 
total column density. Breaks have been found in the distribution of Nhi of absorbers from the 
Lyman-a forest to the DLS (Petitjean et al. 1993; Storrie-Lombardi, Irwin, & McMahon 1996; 
Storrie-Lombardi & Wolfe 2000). We shall show that this behavior of Nhi is compatible with a 
single power law distribution for the total hydrogen column density, once the change in ionization 
fraction is taken into account. In fact it is the change in slope of the Nhi distribution from the 
LLS to DLS region which enables us to derive a value for X. 

A similar phenomenon is found in the outskirts of today's galaxies where a sharp decline of the 
HI column density over a narrow range of the galaxy radius occurs. One can explain this occurrence 
in term of an HI-HII transition zone which takes place when the column density gets sufficiently 
low that the gas becomes optically thin to the extragalactic ionizing flux. For todays galaxies the 
known rotational velocity constrains the dark matter content which, together with the measured 
radial decline of the HI column density at the HI-HII transition zone, can be used to estimate 
the UV ionizing flux at z = 0, otherwise unobservable. Corbelli & Salpeter (1993) have in fact 
modeled the sharp HI decline observed in M33 by considering a total gas distribution compressed 
by the local dark matter (inferred from the observed rotational velocities), and irradiated by the 



-3- 



extragalactic UV and soft X-ray background. They found a best fit to the observed HI radial 
distribution in the outermost disk and to its sharp decline when the intensity of the background 
radiation at the Lyman-edge is 6 x 10 -23 ergs cm~ 2 s _1 sr _1 . For LLS the situation is reversed 
since the metagalactic ionizing flux at high z is better known than the dark matter content of 
absorbers: from X we shall derive the volume gas density n, which is larger than self-gravity alone 
can produce and therefore gives information on the dark matter gravitational potential. 

Unfortunately the difficulty in measuring the residual flux of LLS at large optical depths limits 
the determination of the HI column density of the absorber above 10 18 cm -2 . For this reason 
we are forced to introduce a new technique for treating large uncertainties in the Njji values. In 
the construction of a database for LLS we include also the high column density systems, known as 
Damped absorbers, whose column density is determined from the scatter of the Lyman-a radiation. 
Both the database and the data treatment are described in detail in Bandiera & Corbelli (2000, 
hereafter Paper II). In section 2 of this paper we derive the number density of LLS and DLS from 
our database and in section 3 we describe radiative transfer processes and the gas stratification 
model for the gas in LLS and DLS. The power law distribution for the total hydrogen column 
density and the H fractional ionizations which fit at best the data at 1.75 < z < 3.25 in the Lyman- 
limit and Damped Lyman-a region will be derived in section 4. Having the H ionization fractions 
for LLS, we can then compute the total gas content of LLS and DLS and estimate their dark matter 
potential. Also given in section 4 is an extrapolation of our fit towards absorption systems of even 
lower HI column density. 



2. Number density of Lyman limit systems and Damped Ly-a absorbers 

From the available literature we have collected data on LLS with tll > 0.4 and on DLS with 
rest frame equivalent width W > 5A for z < 4.7. We have 661 QSOs and for each QSO we have the 
redshift path covered with a given sensitivity, the redshift of the intervening absorbers and their 
estimated HI column density, Njji- We disregard absorbers and paths within 5000 km/s from the 
quasar redshift. The sensitivity of a search is expressed as a lower limit to the HI column densities 
searched. To a given direction it may correspond more than one path if different redshift paths were 
searched with different sensitivities. If no estimates of the absorbers HI column density are available 
(if no residual flux is detected short ward of the Lyman break for example) a lower limit to it is 
given while upper limits may be derived from searches of damped Lyman-a absorption lines. The 
HI column density values obtained from Voigt profile fits, when these are available, are included in 
our database. A comparison of the HI column densities derived directly from the equivalent width, 
with those determined from the Voigt profile fit to the same absorption lines shows that the Voigt fit 
procedure gives a systematically higher value of Nhi- This might be due to an underestimate of the 
equivalent width due to an effective absorption of the quasar continuum radiation by intervening 
forest clouds. Using data where Voigt fits are available we derive an average correction for Njji(W), 
the column density derived from W, and we correct the data whenever no Voigt fits are available 
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as follows: 

logN HI = 6.261 + 0.705 logN HI {W), (1) 

Since the maximum dispersion found in the data used for deriving eq. (1) is ±0.4 in logNni 
we use this value for estimating the errors on logNni whenever we apply the above correction. Our 
compilation of data is described in more detail in Paper II and is available upon request from the 
authors. 

The redshift path density w(z) (Lanzetta et al. 1991), defined as the number of lines of sight 
at each redshift searched for LLS with any sensitivity or for DLS with W > 5, is shown in Figure 
1 for all our data. From our data set we extract the LLS sample and the DLS sample. The LLS 
sample consists of all observations which were sensitive to LLS with tll > 1. The DLS sample 
consists of observations sensitive to absorbers with column density logNni > 20.13. We wish to 
examine the number density of absorbers per unit redshift, dN/dz in the LLS and DLS sample. As 
usual we represent the number density as a power law of the form 

N{z) = N (l + Z y (2) 

In order to avoid binning the data in redshift intervals, which can generate errors in the estimate 
of the evolutionary trend, we use the maximum likelihood method described by Storrie-Lombardi 
et al. (1994) to estimate iVo and 7. For the LLS sample the maximum value of the Likelihood is 
found for: 

7 = 1.66 ± 0.3 N = 0.23 (3) 

where the errors quoted are the > 68.3% confidence limits. The > 99.7% confidence limits for 7 
are 0.80 and 2.64. These values are similar to the ones obtained by Stengler-Larrea et al. (1995), 
For the DLS sample we find 

7 = 0.86 ±0.3 iV = 0.18 (4) 

The > 99.7% confidence limits are 0.11 and 1.72. The value of 7 is consistent but slightly lower 
than what has been found by Storrie-Lombardi & Wolfe (2000) for higher column density systems. 
If we restrict our DLS sample to only those systems with have a measured logNni > 20.3 we obtain 
the same values of iVo and 7 as given by Storrie-Lombardi & Wolfe (2000). Figure 2 shows the 
cumulative number of absorbers versus z for both samples. Over plotted is the expected number 
of LLS and DLS from the Maximum Likelihood estimate. Notice that DLS at z > 1.75 seem to 
require a lower value of 7 than what is given by the Maximum Likelihood estimate over the whole 
z range. In fact by considering only redshifts z > 1.75 the match between the data and the best 
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fit cumulative function for DLS improves: in this case 7 = —0.07 with —1.58 and 1.36 as > 99.7% 
confidence limits. 

For the present paper we shall use data only in the intermediate redshift range 1.75 < z < 3.25 
and for this range the slight difference in evolution implied by eq. (3) and (4) is unimportant. 
We use an average 7 value of 1.0 which is consistent both with the number density evolution of 
DLS and of LLS. We exclude from the detections absorbers which were declared as "non damped" 
if they were not detected as LLS or for which no LLS searches have been performed. For "non 
damped" Ly-a absorbers which were detected as LLS the Nhi is set to the value derived from tll 
when this is available; otherwise we shall use the lower limit to tll for a lower limit to Nhi and 
the column density value derived from W as upper limit. 



3. The N H i± - N H ± curve 

Our numerical code solves for the ionization, chemical, and hydrostatic equilibrium of a gaseous 
absorber which receives photons from a background radiation field. It computes the HI column 
density perpendicular to the plane, Nhi±, in a plane parallel geometry when the total hydrogen 
column density perpendicular to the plane is Nh± and the density stratification in the vertical 
direction y is in the following exponential form: 

p{y) =/>oexp(- (5) 

Pq = Mg/(\f2nh) is the central gas density for a given total mass surface density per cm~ 2 of the 
gaseous component, M g . eq. (5) is the exact solution for the Poisson and first moment equation for 
an isothermal gaseous disk in equilibrium in a dark matter halo potential in the limit of a negligible 
self gravity and of negligible changes of the dark matter density with height above the midplane 
(within the gaseous extent, e.g. Spitzer 1942, Maloney 1993). In this case the scale height can be 
written as 



1 _ ^G Pdm _ (2nG) 2 M 2 dm 
h 2 " cl cj 



(6) 



c s is the sound speed (computed using T, the mass averaged temperature of the slab), p^ m and 
Mdm are respectively the dark matter volume density and the dark matter surface density between 
+h and -h in the slab. In the absence of dark matter instead eq. (5) is only an approximation to 
the self gravitating isothermal gas layer solution (e.g. Spitzer 1942) 

9/ V \ 1 2ir 2 G 2 M 2 



sg "sg 
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In this case the gas central density is po = M g / (2\/2h sg ) . If we approximate the above exact 
solution with eq. (5), given the same central volume density po an d surface density M g , we should 
use the following expression for the vertical scale height 



h 2 Ah*,- 4 W 

It will be shown in a subsequent section, that eq. (5), with the scale height h given by eq. 
(8), is a good approximation to the exact self gravitating gas layer solution for the purposes of this 
paper. We can therefore write a general expression for h if we use eq. (5) to describe the vertical 
stratification of the gas when both self gravity and dark matter are taken into account: 



1 16G 2 M g 2 _ ^ 



(-!+'/) -TTT^ (9) 



The compression factor 77 for the volume density is in general considered as the effect of gravity 
from matter other than the gas (dark matter, stars or brown dwarfs which coincide with the gas 
location). We write 77 as 



V = Vo^ (10) 

c s being the value of c s for T = 10 4 K and for a mean gas mass per particle \x = 1.3. In this paper 
we shall consider a set of models, each corresponding to a different 770 value, and in each model 
we keep 770 constant as we very M g while c s is computed from the radiative transfer and energy 
equation. A constant 770 value for 7/0 3> 1 implies 



Pdm oc M 2 g (11) 

The meaning of the above proportionality are not immediate, but for spherical isothermal dark 
matter halos eq. (11) implies a constant ratio between the gaseous and dark total surface densities. 
77 instead is proportional to the the dark matter surface density between +h and —h and since h 
scales linearly with c s , 77 will also scale with c s . 

The equilibrium fractional ionizations of H, He, Hell are found step by step through the slab 
taking into account the diffusion of photons from the recombination of H, He, Hell and the sec- 
ondary electrons produced by the harder photons. The temperature at each step through the slab 
is given by balancing photoionization heating with the cooling from hydrogen and helium gas (col- 
lisional ionizations, recombinations, lines excitation and free-free processes), from metal lines and 
from H2 cooling. After each iteration the mass averaged temperature T is computed consistently 
and is used to determine the gas vertical dispersion and stratification. The ionization states of 
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metals (C,0,Fe) are computed step by step consistently with the photoionization rates and charge 
exchange reactions; very high ionization states (i.e. potential energies E 3> 54 eV) are not consid- 
ered. For metal abundance we use Z = 0.02Z Q throughout this paper but we shall discuss briefly its 
possible variations. H2 fractions, although small, are computed via primordial chemical reactions 
(Corbelli, Galli & Palla 1997). 

From the ionization-recombination balance we know that for a small neutral fraction the ratio 
Nh±/Nhi± depends on the ratio between the ionization and recombination coefficients and varies 
inversely with the gas volume density. We define X as 

X = log^^ for N H i = 1.6 x 10 17 cm" 2 (12) 

Notice that X is defined at the line of sight value of Nhi = 1-6 x 10 17 cm -2 and therefore the 
corresponding Nhi± value depends on the thickness ratio h/R. We derive for ijq > 1 the following 
relation 



X ~ 2.94 + 0.36 log— (13) 
Vo 

where Jl is the intensity of the background flux at 912 A written in units of 10 -22 ergs cm' 2 
s _1 Hz -1 sr _1 . For our redshift range (1.75 < z < 3.25) we use a constant Jl and a flux spectrum 
computed for z = 2.5 by Haardt & Madau (1996, hereafter H-M) for emission by quasars and 
absorption by intervening clouds. A nearly constant ionization rate between z = 2 and 4 is also 
required by comparison of results from hydrodynamical simulations of structure formation with 
the measured opacity of the Lyman-a forest. The H-M flux intensity at z — 2.5, Jl — 5.5, agrees 
sufficiently well with other estimates (Rauch et al. 1997; Giallongo et al. 1996) and we shall use 
Jl = 5.5 throughout. The compression factor r/o is unknown a priori and will be determined by 
trial and error. Hence 770 (or equivalently X), and a are parameters to be determined. 

In Figure 3 we show three curves for Njji± as a function of Njj± for logJi/^o =—1-9, —0.4, 0.7, 
which correspond to X = 2.3, 2.8, and 3.2 respectively. All curves in Figure 3 are characterized by 
3 regions: a first region at high column density where Nh± — Nhi±, a second region where for a 
small change in Njj± the HI column density changes very rapidly, and a third region, towards the 
bottom part of the plot, where the log-/V#j_-log./V//j_L relation is again linear with a slope smaller 
than unity and which in general depends on the variations of r/ with Nh±- For the models in this 
paper, we assume that 770 in eq. (9) does not vary with gas surface density M g (see section 3). The 
factor Jl/(1 + rjo) is what determines the critical column density where the HI — > HII transition 
occurs i.e. the column density where the steep slope of the second region starts. The rapidity of 
the decrease of Njji± in the second region depends also on metal abundances. When Z increases 
above 0.05 Z & metal line cooling becomes quite important in that region and the transition gets 
less sharp as Z increases. The shape of these curves is quite important because it determines the 
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Fig. 3. — The log value of Njji± as a function of Nhj_ for slabs irradiated by the H-M flux at 
z = 2.5 and for logj^/^o = —1.9,-0.4,0.7 (from left to right). The dotted vertical lines indicate 
Njj± = 9.5 x 10 16 cm~ 2 corresponding to Nhi = 1.6 x 10 17 cmT 2 for h/R = 0.2. 
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slope of the distribution function around the Lyman-Limit region. 

4. Prom the observed Nhi distribution function to the total amount of gas 

distributed in LLS and DLS. 

Several papers on the HI column density distribution of absorbers (Tytler 1987; Lanzetta et 
al. 1991; Petitjean et al. 1993) have shown that a power law N H ^ 5 fits the data from 10 13 cm~ 2 to 
10 21 cm" 2 approximately, but not well enough to satisfy statistical tests such as the Kolmogorov- 
Smirnov test. In this paper we shall use the deviations of Nhi distribution from a power law to 
determine X, after assuming that the total gas column density Nh has a power law distribution 
function of the form 

g(N Hi _) = K{l + zVN-l (14) 

In order to derive the Nhi distribution from g(Nj{±) we must orient the absorbers randomly in the 
plane of the sky, assume an axial ratio h/R for the slab, and apply the Nh_l — Nhi± conversion 
factor. For our model fitting we shall use Z = O.O2Z , X independent of redshift and an axial 
ratio of h/R = 0.2, since clouds are likely to be neither spherical nor thin disks. Results essentially 
do not depend on h/R for h/R < 0.5. We derive a and X by comparing the resulting f(Nni) 
for 7 = 1.0 with the data present in our compilation at Njji > 1-6 x 10 17 cm" 2 . K is set by the 
normalization condition for g, based on the observed number of absorbers. 

We emphasize that it is not possible to present the data relative to the distribution function in 
a model independent way, due to LLS with undetermined tll • A deterioration of the available data 
might result from the operation of binning in Nhi m order to render straightforward the comparison 
with the model distributions. Instead of binning the data, we match the model distribution to the 
individual detections. The details of the fitting procedure are given in Paper II; we underline here 
the main characteristics. The procedure is rather similar to that used by Storrie-Lombardi, Irwin, 
& McMahon (1996), but implementing the algorithm in order to take into consideration also the 
uncertainty in the determination of any single value of Nhi- Large observational errors are included 
by leaving undetermined the "real" position of each event. We determine a and X by a Maximum 
Likelihood analysis to the projected HI column density distribution / fixing the real position of 
each event to the measured value of Nhi when this is available; otherwise we use in the Likelihood 
the integral of / between the maximum and minimum value of Nhi- We normalize the theoretical 
distribution such as to give a number of detections with Nhi > 1-6 x 10 17 cm -2 equal to the 
observed one. Two maxima for the Likelihood are found: 



X = 2.82 a = 2.70 K = 1.2 x 10' 



(15) 



-12- 



X = 2.74 



a = 2.57 



if = 1.8 x 10' 



,31 



(16) 



The > 68.3%, > 95.5%, and > 99% confidence levels in the X — a plane are shown in Figure 
4 where the filled dots indicates the location of the maxima as in eq. (15) and (16). The self 
gravitating gas solution (rj = 0, a = 4.63, X = 3.33) lies well outside the > 99% confidence level 
(we would need the 99.999% confidence level to include it) and therefore it is not consistent with 
the data. We have also checked that a similar conclusion holds if we use the exact self gravitating 
gas solution, as given by eq. (7), in deriving the Nhi±_ — Nh± relation. For this case both X and 
the best fit a value are within 1% of the values obtained using eq. (5) and r\ = for the vertical 
gas stratification. 

In Figure 5 we compare the observed value of the cumulative function with the theoretical ones 
derived from the integral of the projected HI column density; we show results for the two best fit 
models (the two maxima in Figure 4) and for two models corresponding to the highest and lowest 
X values on the > 95.5% confidence level of Figure 4. For points which have no defined Nhi, i-e. 
which have large errors, we then compute their best distribution for a given / by spreading the 
data in the allowed range of Nhi according to / weighted with the redshift path. Between all the 
possible permutations of points with undetermined Nhi we then choose those which satisfies best 
the [/-test on the deviations between the observed and the expected cumulative functions over the 
error interval, R, and over Nhi, C (see Paper II for more details and for the use of a numerical 
simulation to proof the validity of this approach). For X — a inside the 99% confidence level of 
Figure 4 the K-S tests on R and C are satisfied to the 99.9% level. 

We have proved that the gas distribution between the LLS and the DLS region follows a single 
power law with index a > 2 if the ionization level is such that less than 1% of the total gas is 
neutral when Nhi = 1.6 x 10 17 cm~ 2 . There is no need of a distribution more complicated than a 
power law once one takes into account ionization effects. Our results on a and X still hold even if 
we do not include in the data set the damped lines with 5 < W < 10 A or if we exclude a certain 
percentage of these lines due to possible blending with smaller lines. 



We shall discuss here some results relative to the best-fit values of X and a as given in eq. 
(15) and in eq. (16), and to the two most extreme values of X on the > 95.5% confidence level, 
namely: X = 3.07 (770 = 2.3), a = 3.32 (hereafter 2cr-up); X = 2.53 (770 = 75), a = 2.27 (hereafter 
2(j— low). The distribution function for the total column density Njj± = Nhj_/W 20 cm -2 can be 
written as 



4.1. The total gas content of LLS and DLS 



g(N H ±) = K(l + z)i\T 



H± 



—a 



(17) 
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Fig. 4.— The > 68.3%, > 95.5%, and > 99% confidence levels in the X - a plane. 
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Fig. 5. — The theoretical cumulative distribution functions (continuous lines) for our two best fit 
models ((a) a = 2.70, X = 2.82 and (b) a = 2.57, X = 2.74) and for two points over the > 95.5% 
confidence level ((c) a = 3.32, X = 3.07 and (d) a = 2.27, X = 2.53). For comparison the observed 
HI cumulative functions (dotted lines) are shown. 
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In Figure 6, in arbitrary scale, the continuous lines show log f(Nni), the HI distribution 
function, for the best-fit values of X and a as given in eq. (15) and (16), and for the 2a— low and 
the 2cr-up models. Our data for Nhi > 1.6 x 10 17 cm' 2 is in five large bins just for the purpose of 
presentation. 

Nh± x q(Nh±) can be integrated over a range of Nh±, say between Nn±,i and Nn±, n , to 
estimate the mass density of hydrogen atoms in gas clouds with an average HI column density 
along the line of sight between Nhi,i and Nhi, u - The comoving cosmological H+He gas density at 
z = 2.5 can be written as: 

IWOto - 2 *Z7ll% ^ - "»tn) (18) 

5 depends on the cosmological model and is a function of z,QmMa- For a standard Friedmann 
Universe in which qo = 0, S = 1 and we shall use this value for the rest of this section (for VLm = 0.3 
and Qa = 0.7 instead S depends on z and is close to zero at z ~ 2.5). 

In Table 1 we give the values of f2 ffas (i)/i6o, for i = 1,2,3 which is the total gas density in 
the Universe at z ~ 2.5 due to absorbing clouds whose HI column density projected along the 
line of sight is between Nhi,i and 10 22 cm -2 . We shall consider Nhi,i = 10 14 (see section 4.2), 
Nhi, 2 = 1-6 x 10 and Nhi,3 = 1-3 x 10 20 cm" 2 . For each corresponding value of Nn±i we give 
Xi, the log ratio of total to neutral gas column density. Results are given both for the best fitting 
models and for the two most extreme values of X on the > 95.5% confidence level of Figure 4. In 
the Table we also show the gas scale heights for Nhi,2 = 1-6 x 10 17 and Nhi,3 = 1.3 x 10 20 cm" 2 , 
and values of ^dark(i)- 

For Qdark in the regions coinciding with the gas the factor ( to be substituted into Table 1 is 
C ~ 1, independent of assumptions on cloud size and relative distribution of dark matter and gas. 
For rotating disks embedded in spherical dark halos one can compute the contribution of the total 
dark matter surface density to the cosmological matter density. This contribution associated with 
DLS or LLS systems depends on the rotational velocity V and is given by &dark using £ rj V/c s . 
This factor may be close to one for dwarf clouds, but £ S> 1 would hold for giant disk proto- galaxies. 
For the range of uncertainties in Table 1, the value of Q gas (2), the total contribution of LLS plus 
DLS, varies little but 770 has a large spread. Consequently, Qdark and the gas scale height h also 
have a large spread. Note that h(3) is particularly small for the 2cr-low limit. 



4.2. Conjectures on the Lyman-a forest 

The Lyman-a forest clouds may not be directly related to LLS, but we can explore the con- 
sequences of two assumption, namely that eq. (17) extends into the Lyman-a forest region with a 
compression factor 770 in eq. (9) that has the same constant value as for LLS and DLS. The filled 
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Fig. 6. — HI column density distribution functions between z = 1.75 and z = 3.25 for X = 
3.07,2.82,2.74,2.53 and a = 3.32,2.70,2.57,2.27 respectively (from top to bottom). Bins with 
open circles have data from our collection while filled triangles are data from the compilation of 
Petitjean et al. (1993). Scale is arbitrary for each theoretical model but data points have been 
scaled accordingly. The long dashed line indicates the extrapolation of the distribution function 
towards the forest region using the same r/o value as for the LLS and DLS. 



-17- 



triangles in Figure 6 are the observed data for the Lyman-a forest from Petitjean et al. (1993) 
and Hu et al. (1995). The dashed curves show the extrapolations of the distribution functions 
keeping a and rjo constant. Note that these extrapolations, without any additional parameters, 
give surprisingly good fits, especially the best fit model given by eq. (15) and the 2<r-up model (the 
2cr— low model and the best fit model given by eq. (16) favour a different mechanism of confinement 
in place for the very low column density clouds, like external pressure, which would steepen the 
distribution by keeping the fractional ionization of H independent from the cloud column density) . 

If Lyman-a forest clouds were physically different from LLS and DLS, one might have expected 
a change in a or rfo and hence a much worse fit of the dashed curves in Figure 6. The good fit 
makes it even more likely that a and rjo are constant from LLS to DLS, which represents a much 
smaller range in N h- 

5. Summary and Discussion 

The results of the model fit to LLS and DLS presented in this paper underline new aspects in 
the column density distribution of Lyman-a absorbers: 

• The excess of systems in the damped region, suggested by the Nhi data, is naturally explained 
in terms of a sharp transition from a highly ionized to a highly neutral gas distribution and it does 
not requires any break in the distribution of the total gas column density. 

• The total gas column density distributions, g(Njj±), which best fits the data for 1.75 < z < 
3.25 in the LLS and DLS region can be described by a power law of index —a with 2 < a < 3.7. 
This has the important consequence that low column density systems contains more mass than 
high column density systems. 

• We have tested that our results, which are relative to the redshift bin 1.75 < z < 3.25, 
depend weakly on the data selection and on the number redshift evolution. They also do not show 
a strong dependence on the assumed thickness to diameter ratio of the slab or from its metallicity 
if Z < 0.05Z o . 

• The gas fractional ionizations increase with decreasing column density and data are best fitted 
when hydrogen fractional ionizations are of order ~ 0.002 at tll = 1 i-e. X ~ 2.8. Gas fractional 
ionizations for absorbers in a background radiation field depends on the gas volume densities and 
temperature. The model presented in detail in this paper considers the gas self gravity and a dark 
matter potential in which the dark matter surface density inside one gas scale height is proportional 
to the gas surface density (constant r]o). However the range of a and X values compatible with the 
data are similar if one considers instead a model where 770 varies with Njj (e.g. if rjo x Nh is kept 
constant instead). Table 1 summarizes our main results; the values of 770 for more general models 
refer to LLS column densities which are no longer opaque to the UV ionizing radiation. 

Although our Likelihood analysis allows power law indices a in a rather wide range, 2 to 3.7, 
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other arguments make the upper half of this range more likely: a > 2.7, r/o < 13, X > 2.8. For 
larger values of r/o the gas scale height h(3) for DLS is uncomfortably small and the LLS contribution 
^cZarfc(2) to cosmological matter density is uncomfortably large (see Table 1, the £1 values given are 
for qo = 0; values for Om ~ 0.3 and Q\ = 0.7 are larger by a factor ~ 2.5). For the total dark 
matter mass, the factor ( for Sldark is of order V/W km s _1 , which suggests V < 100 km s _1 on 
average for LLS and DLS at redshifts z ~ 2.5. This finding is in agreement with the rotational 
velocities of absorbers halos predicted by Valageas, Schaeffer & Silk (1999) and with models that 
propose that a large fraction of Lyman-a absorption systems originate from small low luminosity 
systems (Abel & Mo 1998; Rao & Turnshek 2000; Haehnelt, Steinmetz & Rauch 2000). However 
Steidel, Dickinson & Persson (1994) have shown that small rotational velocities are not required by 
low redshifts Mgll absorption systems with W > 0.3 A. These have the same average density per 
unit path of redshift as LLS and seem to be associated with normal galaxies whose halo extends for 
~ 40 kpc. This opens the issue about possible evolutionary scenarios for LLS, which can be fully 
addressed only once we know the redshift evolution of galactic halos and of the galaxy luminosity 
function. The possible role of faint galaxies as those in pairs close to bright galaxies (Churchill 
2000) should also be taken into account. Futhermore it would be extremely useful to have a larger 
statistical sample of LLS and DLS at z < 1.75 and at z > 3.25 which, together with a better 
knowledge of the evolution of the background ionizing radiation field, allows a determination of a 
and i] at these redshifts (see Paper II for an attempt with the actual data set). 

• Although Lyman-a forest clouds, with their much smaller column density, might be physically 
different and have different values of a and r/o, extrapolations using constant a and 770 values as for 
LLS and DLS give surprisingly good fits. 

We cannot make definitive statements about the low density Lyman-a forest clouds because 
of possible dynamic deviations from hydrostatic equilibrium and possible pressure confinement. 
Nevertheless, the fact that a is appreciably larger than 2 near LLS for z ~ 2.5 suggests that forest 
clouds contained more gas than LLS, which in turn contained more gas than DLS. Several papers 
have shown that the gas content in DLS at z ~ 2.5 is below the mass content in galaxies in the 
local Universe (Storrie-Lombardi & Wolfe 2000), Q s+g h^o ~ 3 — 7 x 10~ 3 and that the stellar 
mass density in stellar systems at z ~ 2.5 is very small (Madau, Pozzetti, Dickinson 1998). It is 
therefore likely that gas in forest clouds and LLS at z ~ 2.5 has collapsed and contributed to stars 
in present day galaxies. 

We are grateful to P. Madau for providing us numerical results on the background radiation 
field, to the referee and to Dr. Bothun for very useful comments to the original manuscript, and to 
D. Chernoff for helping with the Arcetri-Cornell connection. One of us (RB) acknowledges partial 
support by NSF-PHY94-07194. 
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Table 1. Ionization and gas content of Lyman-a absorption clouds. 
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